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Abstract 

Self-similar solutions of a collapsing perfect fluid and a massless scalar field with kinematic self- 
similarity of the first kind in 2+1 dimensions are obtained. Their local and global properties of the 
solutions are studied. It is found that some of them represent gravitational collapse, in which black 
holes are always formed, and some may be interpreted as representing cosmological models. 

1 Introduction 

One of the most interesting problems in gravitation theory is the study of the relation that exists between 
the critical phenomena and the process of black hole formation. The studies of non-linearity of the 
Einstein field equations near the threshold of black hole formation reveal very rich phenomena [l]-[3], 
which are quite similar to critical phenomena in Statistical Mechanics and Quantum Field Theory [4]- [5]. 
In particular, by numerically studying the gravitational collapse of a massless scalar field in 3 + 1- 
dimensional spherically symmetric spacetimes, Choptuik found that the mass of such formed black holes 
takes a scaling form, Mbh = C{p) (p — p*) 1 , where C(p) is a constant and depends on the initial data, 
and p parameterizes a family of initial data in such a way that when p > p* black holes are formed, and 
when p < p* no black holes are formed. It was shown that, in contrast to C(p), the exponent 7 is universal 
to all the families of initial data studied. Numerically it was determined as 7 ~ 0.37. The solution with 
p = p* , usually called the critical solution, is found also universal. Moreover, for the massless scalar field 
it is periodic, too. Universality of the critical solution and exponent, as well as the power-law scaling 
of the black hole mass all have given rise to the name Critical Phenomena in Gravitational Collapse. 
Choptuik's studies were soon generalized to other matter fields [6, 7], and now the following seems clear: 
(a) There are two types of critical collapse, depending on whether the black hole mass takes the scaling 
form {Mbh) or not. When it takes the scaling form, the corresponding collapse is called Type II collapse, 
and when it does not it is called Type I collapse. In the type II collapse, all the critical solutions found 
so far have either discrete self-similarity (DSS) or homothetic self-similarity (HSS), depending on the 
matter fields. In the type I collapse, the critical solutions have neither DSS nor HSS. For certain matter 
fields, these two types of collapse can co-exist, (b) For Type II collapse, the corresponding exponent is 
universal only with respect to certain matter fields. Usually, different matter fields have different critical 
solutions and, in the sequel, different exponents. But for a given matter field the critical solution and the 
exponent are universal. So far, the studies have been mainly restricted to spherically symmetric case and 
it is not clear whether or not the critical solution and exponent are universal with respect to different 
symmetries of the spacetimes [8, 9, 10]. (c) A critical solution for both of the two types has one and 
only one unstable mode. This now is considered as one of the main criteria for a solution to be critical, 
(d) The universality of the exponent is closely related to the last property. In fact, using dimensional 
analysis [11]- [14] one can show that 7 = jjy, where k denotes the unstable mode. 
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From the above, one can see that to study (Type II) critical collapse, one may first find some 
particular solutions by imposing certain symmetries, such as, DSS or HSS. Usually this considerably 
simplifies the problem. For example, in the spherically symmetric case, by imposing HSS symmetry the 
Einstein field equations can be reduced from PDE's to ODE's. Once the particular solutions are known, 
one can study their linear perturbations and find out the spectrum of the corresponding eigenmodes. If 
a solution has one and only one unstable mode, by definition we may consider it as a critical solution 
(See also the discussions given in [15]). The studies of critical collapse have been mainly numerical so 
far, and analytical ones are still highly hindered by the complexity of the problem, even after imposing 
some symmetries. 

Lately, Pretorius and Choptuik (PC) [16] studied gravitational collapse of a massless scalar field in 
an anti-de Sitter background in 2 + 1-dimensional spacetimes with circular symmetry, and found that 
the collapse exhibits critical phenomena and the mass of such formed black holes takes the scaling form 
(Mbh) with 7 = 1.2 ± 0.02, which is different from that of the corresponding 3 + 1-dimcnsional case. 
In addition, the critical solution is also different, and, instead of having DSS, now has HSS. The above 
results were confirmed by independent numerical studies [17]. However, the exponent obtained by Husain 
and Olivier (HO), 7 ~ 0.81, is quite different from the one obtained by PC. It is not clear whether the 
difference is due to numerical errors or to some unknown physics. 

After the above numerical work, analytical studies of the same problem soon followed up [18, 19, 20, 
21, 22, 23]. In particular, Garfinkle found a class, say, S[n], of exact solutions to the Einstein-massless- 
scalar field equations and showed that in the strong field regime the n = 4 solution fits very well with 
the numerical critical solution found by PC. Lately, Garfinkle and Gundlach (GG) studied their linear 
perturbations and found that only the solution with n — 2 has one unstable mode, while the one with 
n = 4 has three [22]. According to 7 = pr, the corresponding exponent is given by 7 = l/|fc| = 4/3. 
Independently, Hirschmann, Wang & Wu (HWW) systematically studied the problem, and found that 
the n = 4 solution indeed has only one unstable mode [23]. This difference actually comes from the 
use of different boundary conditions. As a matter of fact, in addition to the ones imposed by GG [22], 
HWW further required that no matter field should come out of the already formed black holes. This 
additional condition seems physically quite reasonable and has been widely used in the studies of black 
hole perturbations [24]. However, now the corresponding exponent is given by 7 = l/|fc| = 4, which is 
significantly different from the numerical ones. So far, no explanations about these differences have been 
worked out, yet. 

Self-similarity is usually divided into two classes, one is the discrete self-similarity mentioned above, 
and the other is the so-called kinematic self-similarity (KSS) [25, 26], and sometimes it is also called 
continuous self-similarity (CSS). KSS or CSS is further classified into three different kinds, the zeroth, 
first and second. The kinematic self-similarity of the first kind is also called homothetic self-similarity, 
first introduced to General Relativity by Cahill and Taub in 1971 [27]. In Statistical Mechanics, critical 
solutions with KSS of the second kind seem more generic than those of the first kind [4]- [5]. However, 
critical solutions with KSS of the second kind have not been found so far in gravitational collapse, and 
it would be very interesting to look for such solutions. 

Besides, recently we have studied gravitational collapse of perfect fluid with kinematic self-similarities 
in four-dimensional spacetimes [28], a subject that has been recently studied intensively (for example, see 
[29]-[32] and references therein.). In this paper, we shall investigate the same problem but in 2+1 gravity 
[33, 34] and including a massless scalar field [35] [36]. The main motivation of such a study comes from 
recent investigation of critical collapse of a scalar field in 3D gravity [16, 17, 18, 19, 20, 21, 22, 23]. 

We shall present in this work the self-similar solutions of a collapsing perfect fluid and a massless 
scalar field with kinematic self-similarity of the first kind in 2+1 dimensions. Their local and global 
properties of the solutions are studied. Specifically, in section II we present the self-similar Einstein's 
field equations, and in section III we show the general solutions obtained. In section IV we present the 
global properties of the solutions and in Section V we conclude our work. 

2 Field Equations 

The 2+1-dimcnsional spacetimes with circular symmetry and zero-rotation are described by the metric 

ds 2 = j ab (t, r)dx a dx b + g ee (t, r)d9 2 , (1) 
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where {x a } — {t, r}, (a, b — 0, 1), and 8 denotes the angular coordinate, with the hypersurfaces 8 = 0, 2ir 
being identified. The metric is invariant under the following coordinate transformation, 

t = t(t>,r>), r = r(t',r'). (2) 
On the other hand, the energy-momentum tensor for a perfect fluid is given by 

T^{p,P) = {p + p)u^u v - pg^ , (3) 

where is the velocity of the fluid, p and p are its energy density and pressure. The energy-momentum 
tensor for a massless scalar field is given by 

r F W = 4> tli <t>, v - \gp.v<i>., a ^ a , (4) 

where ( ) jQ , = d( )/d a . We choose the coordinates such that 

n-ii = Vdoo <5° , and g 01 = , (5) 

which implies that the coordinates are comoving with the perfect fluid. Then, the metric (1) can be cast 
in the form 

ds 2 = l 2 ^ tr Ut 2 - e 2 *('^rfr 2 - r 2 S 2 {t, r)^ 2 | , (6) 

where I is an unit constant with the dimension of length, so that all the coordinates {x a } = {t, r, 8} are 
dimensionless. We construct a projector operator h^ v by 

= g^v - u^Uy , (7) 

from which we find h a ^u a Ufj = 0. Once defined the projection operator h^, following Carter and 
Henriksen [25, 26], we define the kinematic self-similarity by 

= 2h tlv , = -au^ , (8) 

where denotes the Lie differentiation along the vector £ M , a is a dimensionless constant. When a = 0, 
the self-similarity is said to be of the zeroth kind, when a = 1 it is of the first kind (or homothetic), and 
otherwise it is of the second kind. In this work we shall consider only solutions with self-similarity of the 
first kind (a = 1). Applying the above definition to the metric (6), we find that 

r) = $(x) , r) = , S(t, r) = S(x) , (9) 

where the self-similar variable x and £ M are given by 

x = \n(^j, (10) 

and 

Before looking for solutions of the Einstein field equations, we must take into account that for the 
metric to represent circular symmetry some physical and geometrical conditions must be imposed [37]- [46]. 
We impose the following regularity conditions for the gravitational collapse: 

(i) There must exist a symmetry axis, which can be expressed as 

H=\^w9^\ 1/2 ^0, (12) 

as r — > 0, where we have chosen the radial coordinate such that the axis is located at r — 0, and £^ is 
the Killing vector with closed orbit given by £?g\d a — 8$. 

(ii) The spacetime near the symmetry axis is locally flat, which can be written as [47] 

K, a 1l iP g a P -> -1 , (13) 
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as r —* 0. Note that solutions failing to satisfy this condition sometimes are also acceptable. For example, 
when the right-hand side of the above equation approaches a finite constant, the singularity at r = 
may be related to a point-like particle [48]. However, since here we are mainly interested in gravitational 
collapse, in this paper we shall assume that this condition holds strictly at the beginning of the collapse, 
so that we can be sure that the singularity to be formed later on the axis is due to the collapse. 

(iii) No closed timelike curves. In spacetimes with circular symmetry, closed timelike curves can be 
easily introduced. To ensure their absence, we assume that the condition 

holds in the whole spacetimc. 

In addition to these conditions, it is usually also required that the spacetime be asymptotically flat 
in the radial direction. However, since we consider solutions with self-similarity, this condition cannot 
be satisfied, unless we restrict their validity only up to a maximal radius, say, r = ro(t), and then join 
them with others in the region r > r (t), which are asymptotically flat as r — > oo. In this paper, we shall 
not consider such a possibility, and simply assume that the self-similar solutions are valid in the whole 
spacetime. 



3 Self-Similar Solutions of the First Kind 

In this section, we study solutions with self-similarity of the first kind. To obtain the desired equations 
we substitute equations (3), (4) and (5) into the Einstein field equations 

G^ = k[T„ v { P , P )+T^(cP)], (15) 

and also using the Klein- Gordon equation 

□0 = 0, (16) 

where n is the Einstein coupling constant, □ = g a(3 V a Vi3 with V Q being the covariant derivative. When 
a = 1, according to equation (10) terms with powers of r can be substituted by the same powers of 
(—t) since r = c x (—t). Then it can be shown that the Einstein field equations in this case become (see 
Appendix A) 



e 



JM) [»*, I -^ + WVe-»[ tf , I + (l + y)( !/ -*, I ) + J(W J ] = ?Po, (17) 



-e 2 (*-*>[^ + y (y-<I>, x + l) + 5(^ = Ipo, (18) 



-e 2 <*-*> [*, xa + - $ , + 1) + fa x + 0, T ) 2 ] + 



e" 2 *[<I>, xx + <M*x-*x-l) + ^,*) 2 ] = -p , (19) 

y,* + (i + y)(tf-*,x)-y$,* + K<M£x + £T) = °. ( 2 °) 

e 2( *- $) [(</>,, + 4>.r), x - - * x - V - l)(0,x + 0,r)] - 

e- 2 *[<A,xx + ($,*-*x + yMx] = 0. (21) 

In order to determinate the metric coefficients and the x dependence of the scalar field completely, 
and to make possible the calculations, we assume the equation of state for the perfect fluid given by 



P = up , (22) 

for any u. 

Wc use for the massless scalar field the general form <p(r, x) = ar + F(x), where a=constant. 

To solve the above equations in the general case is not an easy task. In the following, wc shall consider 
several particular cases. We first show some particular solutions for a dust fluid and later for the general 
case with w / 0. Solutions corresponding to p = p = 0, that have already been considered in reference 
[23], are also shown in section 3.3. 
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3.1 Dust Fluid (p = 0) 

From geodesic equations we find that $(x) = $o, constant. Then, the Einstein equations are given by 





-f(^x + ^) 2 ]-e- 2 *[j/,x + (H 


■2/)(y - 


*x) + f(^,x) 2 ] 


K 

= 2 P0 ' 


(23) 




-e 2(3; -* o) [2/,x + 2/(j/+l) + f(^ 




]-e- 2 *f(^) 2 


= o, 


(24) 


_ e 2(- 


"* o) [*,** + *,*(*,* + l) + |(0 >it 




]+e- 2 *f(^) 2 


= o, 


(25) 




y,x + + 


■*,*) + 


« 4>,x($,x + 4>,t) 


= o, 


(26) 


- $0) [(^x + ^) ; 


* + (*,* + y+l)(£* + 0,r)] -e" 


2 *[^x 


- {^, x -y)4>, x ] 


= . 


(27) 



In the following, we consider several particular cases. 
3.1.1 Solution 1 



*(a:) 


= *0 




= -x + * 


y{x) 


= -1 


S(x) 


= S e~ x 


4>{x,t) 


= 4>o 




e -2*„ 


p 




p 


= 



where $o, So and <f>o are arbitrary integration constants. 
3.1.2 Solution 2 



3.1.3 Solution 3 



*(a:) 


= $0 


*(x) 


[it 


= - a ]/2 X - 


y{x) 




= a \jl- 1 


S{x) 


= S e y x 


4>(x,t) 


= CIT + 00 


p 


= 


V 


= o, 




I 


a 









(28) 



$(#) = $ 

*(*) = ln[cosh(^)]-| + * 

y(x) = -1 

S(x) = S e~ x (29) 

<j>(x,r) = 0o 

g-2*o X — Xq 

P = 2«^(-t)» [l-tanh(— 3— )] 

p = 



(30) 
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3.2 p = up 

Now, let us consider particular solutions for w =^ 0. Thus, we will use the general equations (17)-(21). 



3.2.1 Solution 4 



$(x) 

y{x) 

S(x) 
4>{x,t) 

p 
p 



UJ+1 

-1 
w+1 

S e BI 

00 



+ *0 



-2*„ 



K Z 2 (w + l) 2 (-t) 2 
UJp 



(31) 



3.2.2 Solution 5 



$(x) 

0(x,t) 

P 
P 



= $ 



o 
1 

~2 

Soe- l/2 
a t + 4>o 

e -2*o 

4k J 2 (-i) 2 
P 



(1 - 2ko 2 ) 



(32) 



3.2.3 Solution 6 



<E>(x) 

2/0*0 

S(x) 



= ±$ 1 ln|c 2: + 2 * - 1 1 + k a 2 x + $ 

= *(«)-| + *o 

= -1/2 

= a(r - x) + 0o 



P 



e- 2 *°(2«; a 2 -l±2$i) 



(33) 



4 K ;2(_£)2(1-Ka 2 ) r 2rea 2 | e x+2* _ ]_ | ±2*i 



with + if x + 2*o > and - if x + 2* < 
3.2.4 Solution 7 



$0*0 
*(x) 

2/0*0 

S(x) 
<x,t) 



$0 

-x + * , *o = $o 



P = 



1 + Kfl 2 

o(r + yx) + 4> 

e " 2 " r (i-y)- 1 + 2y 



(34) 



«Z 2 (-/) 2 



2y 



6 



p 



l + y 

1-2/ 



± 



l + y 



k y 



1 < y < 



3.3 p = p = 
3.3.1 Solution 8 



3.3.2 Solution 9 



3.3.3 Solution 10 



3.3.4 Solution 11 



2/(a;) 
0(z,t) 



*(a:) 
0(x,r) 



$ 

I[tanh(^)-1] 
5oe-^ 2 cosh(^-^) 



$(z) 
2/0*0 



$0 

x 

-2+*° 
-1/2 

p = 
1 



$(x) 

2/0*0 

S(x) 
4>{x,t) 

p 



= ± 



$0 

-X + * , *0 = $0 

-1/2 
See'*/ 2 

o(t - z/2) + 0o 

p = 
1 



1 - 2ko 2 



(35) 



(36) 



(37) 



ln|e a:+2 * - 1| + na 2 (x + 2* ) + $o 



$ (a; ) - - + * 

-1/2 

5oc-/ 2 

o(r - x) + 0o 

p = 



(38) 



3.3.5 Solution 12 



$(x) 


x 
2 




X 

= -- + *„, * 
z 


y{x) 


= -1/2 


S(x) 


= S e-*/ 2 


4>{x,t) 


= o(t - z/2) + </> 


p 


= p = 


a 







(39) 



4 Properties of the Self-Similar Solutions 

In this section we shall analyze the global structure for the solutions found in the last section. 
4.1 Dust Fluid 

Let us first consider the solutions for p = 0. 
4.1.1 Solution 1 

From equation (28), we find that the corresponding metric can be written as 

{-tf 



ds 2 = l 2 l e 2 *°dt 2 



The geometrical radius is given by 



K = lS re- x =lS (-t) 



(40) 



(41) 



Without loss of generality, we assume that So > 0. Since 1Z is a function of t only, it is easy to see that 
IZj is always positive, since —t > 0. The whole spacetime is trapped, as one can see from the outgoing 
and ingoing null geodesies expansions, which are now given by the equations (see Appendix A) 



-*o 



and 9„ 



(42) 



2g(-t) " 2/(-t) ' 

from which we have that 9i6 n > . 

From equation (28) we see that the space-time is always singular when (— t) — ► 0. On the other hand 
the following expression, 



R, a R,i3 9 



a/3 



l 2 S 2 c- 2x 



-2*n 



(-ty 



- 



l 2 S 2 c- 2 ^ 



(43) 



is always nonzero. In this solution we do not have an apparent horizon but a singularity at t = 0. Thus, 
it may be interpreted as representing a cosmological model. 



4.1.2 Solution 2 

In the case of equation (29), the metric reads 

p 2*o+x 

ds* = i l { e 2 ^dt 2 - —. 



(-*) 



dr 2 - S 2 (-t) 2 d6 2 



(44) 



The outgoing and ingoing null geodesies expansions 
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Oi 



and 6 n 



(45) 



2g(-t) 2f(-t) 

are such that &i9 n > . Thus, there is no apparent horizon. 
The geometrical radius is given by 

R = lS re- x =lS (-t) , (46) 
where it is assumed that So > 0. On the other hand, we see that the right hand side of the expression 

„2 



R a R P g aP = i 2 s 2 c- 2 * 



-2*0 



H) 2 







; 2 So 2 e- 2 *° 



(47) 



is always nonzero. 

This case is similar to that of Solution 1. Thus, this solution, as in Solution 1, may be interpreted as 
representing a cosmological model. 

4.1.3 Solution 3 

We can see that this solution is equal to Solution 9, which is a solution studied by Hirschmann, Wang & 
Wu [23]. 

4.2 p = up 

Now, we consider the solutions of equations (31)-(34) for ic^O. 

4.2.1 Solution 4 

From equation (31) we see that for < ui < 1 the space-time is the same as that obtained by Miguelote 
et al. [34], but for a constant scalar field which, without loss of generality, can be 4>o = 0. 

4.2.2 Solution 5 

From equation (32), we find that the corresponding metric can be written in the form 



<ls 2 --- lH V- h .I! 2 - ( -Jl^o dr 2 + r 2 5 2^2 } 



The geometrical radius is given by 



H = lrS e- x/2 = lS r 1/2 (-t) 1 



/2 



(48) 



(49) 



We assume that Sq > 0, since 1Z > 0. 

Let us now investigate the behavior of the ingoing and outgoing expansions of the null geodesies. 
From equations (32), (70) and (71) we find that 



and 



where 



0i = 



4ff(-t) 



1/2 



On — — 



4/(-t) 



1/2 



> if r < v ah 
= if r = r ah 
< if r > r AH 



< , for any r G [0, oo) 



(50) 



TAH - C 2 (*°-*o)(_ t) 



(51) 
(52) 
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is the apparent horizon. From the above equations we see that 

> if r > tah — trapped 

= if r = tah — marginally trapped (53) 
< if r < tah — untrapped 

In order to satisfy the regularity condition (13) we find that 

So = 2e*° . (54) 

Thus, we conclude that we have a black hole. 

4.2.3 Solution 6 

Due to the discontinuities in p and 9 n at the horizon, this solution does not represent a physical system. 

4.2.4 Solution 7 

From equation (34) we find that the metric is given by 

ds 2 = Z 2 {e 2 *° dt 2 - ^e 2 * dr 2 - S 2 ^-^ dfl 2 j . (55) 
The geometrical radius is given by 

r i+y 

U = lrS ^ = ISo-— , (56) 

where we can see that 1Z — > when r — > 0, since — 1 < y < because a is real. We also assume So > 0. 
The outgoing and ingoing null geodesies expansions are given by 



/3w- 1 



, , ■ ' «' ;i w>1 /3 (?/>-l/2) 

9t = — — =0 if u = 1/3 (y = -1/2) (57) 

2gr V^ + iy ( <0 if w<1/3 (y < -1/2) 



and 

e -*o«o 



< for any r £ [0, oo) . (58) 



2/r 

Thus, this solution does not present an apparent horizon and it represents a cosmological solution. 

4.2.5 Solutions 8 

As can be seen in equation (35), this solution, corresponds to a flat space-time. 

4.2.6 Solutions 9, 10, 11 e 12 

Since these solutions have already been studied by Hirschmann, Wang & Wu [23], we do not present any 
analysis of them. 

5 Conclusions 

In this paper we have obtained self-similar solutions of the Einstein field equations for a collapsing massless 
scalar field and perfect fluid with kinematic self-similarity of the first kind in 2+1 dimensions. The local 
and global properties of the solutions are studied. Since in the general case it is not an easy task to solve 
the system of field equations, we have considered some particular cases. We have found 12 solutions. 
One of them does not represent a physical system (Solution 6) and another one of them represents a 
Minkowski spacetime (Solution 8). Six of them represent solutions already studied by others authors 
(Solutions 3, 4, 9, 10, 11 and 12). Thus, we have obtained three new cosmological solutions (Solutions 1, 
2 and 7) and a new black hole solution (Solution 5). 
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Appendix A - The field equations 

The non-zero components of the Einstein tensor are given by 

Gtt = + + ( 59 ) 

_ e 2(*-$) x 

Grr = (_ t ) 2 + y(y - + 1)] + ^(i + i/)$,x (60) 



Gee = S 2 r 2 \ + $,* + !)] 



+ + $,x($* -**-!)] j (61) 

G tr - ^j^: [i/,x + (i + y)(y -*,»)- »*,x] , (62) 



where 



!/4 = ( ln5 )-- > ( 63 ) 
and the Klein- Gordon equation can be written as 

72 lV,xx ~ ($ x - * x - y - l)(^,x + <f>,r)] 



e -2* 



(-tf 

e -2* 

[<^,xx + (*, x - *,x + y)4>,x] = o ■ (64) 



The non-zero components of the energy-momentum tensor for a massless scalar field and a perfect fluid 
are given by 

Tu = ^(^+W 2 + ^(W 2 + Pe 2$ (65) 



e 2(*-*) 1 

Tee = r 2 S 2 



Trr = ^T^2 (0,x+0,r) 2 + — fcf+pC 2 ' (66) 

-2* 



+ 0.r) 2 -^-5-(^x) 2 +P 



2 {-t) 2 ^ ^' TJ 2 r 



(67) 



T tr = T^y^xW.x + ^.r)- (68) 

In order to put the field equations in a more suitable form for calculations we write the energy density 
and pressure as 



Po{x) , Po(x) 

p = Y^ P =T^- (69) 



For the self-similar solutions of the first kind the outgoing and ingoing expansion of the null geodesies 
and 1Z iCt lZ^g al3 are, respectively, given by [34] 
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<? 1 =2^{(l + y)e-* (x) +ye x -* (l) } , (70) 

On = | (1 + y)e-*(*> - ye 1 *-*™ } , (71) 

where / and 5 are assumed to be positive functions, / > and g > (for a more detailed discussion, see 
appendix B of reference [34]), and 

fc,„W = 5 2 (x)| y 2 ^^ - (l + ») a e- a *W| . (72) 
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